On shape and electrostatics: competing anisotropics in charged colloidal platelets 
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Charged platelet suspensions, such as swelling clays [2-0] 5 disc-like mineral crystallites 0, Hj] or 
£v^| ■ exfoliated nanosheets [fiLLD], are ubiquitous in nature. Their puzzling phase behaviours are nevertheless 
still poorly understood: while Laponite and Bentonite clay suspensions p], [1| form arrested states 
at low densities, others, like Beidellite and Gibbsite 0, 0], exhibit an equilibrium isotropic- nematic 
H-H ' transition at moderate densities. These observations raise fundamental questions about the influence 
of electrostatic interactions on the isotropic-nematic transition and more generally on the organisation 
of charged platelets. Here, we investigate the competition between anisotropic excluded-volume and 
electrostatic interactions in suspensions of thin charged disks, by means of Monte-Carlo simulations. 
We show that the original intrinsic anisotropy of the electrostatic potential between charged platelets, 
obtained within the non-linear Poisson-Boltzmann formalism, not only captures the generic features 
of the complex phase diagram of charged colloidal platelets, but also predicts the existence of novel 
structures and arrested states upon varying density and ionic strength. 

Colloidal suspensions of spherical particles have played a historical role as model systems to elucidate a variety of 
fundamental phenomena such as glass transition and crystal nucleation [H, However, many naturally occurring 
colloids such as tobacco mosaic viruses, DNA strands, swelling clays, and other mineral crystallites are anisotropic. 
Recent breakthroughs in the synthesis of particles with precisely controlled shapes and interactions have sparked re- 
newed interest in understanding the influence of particle anisotropy on phase behaviour [l(J [H| , which is required to 
achieve well-controlled self-assembly and design of novel structures. In the simplest case of axially symmetric particles, 
1 anisotropic excluded-volume interactions lead to formation of liquid-crystalline phases [12] while more complex shapes 
■ result in remarkable structural diversity [11]. Directional interactions between particles, induced by their geometry and 
1 chemical structure, can also lead to anisotropic complex self- assembled structures. When both particle geometry and 
interactions are anisotropic, the interplay, and even competition between the excluded-volume and other (e.g. electro- 
ns! ' static) interactions may lead to the emergence of novel structures and geometrical frustration. Charged platelet and 
nanosheet suspensions, that have been widely studied experimentally and that exhibit an intriguing phase behaviour, 
provide such examples and are the subject matter of this report. For hard-core repulsions only, upon increasing 
platelet density, a sequence of isotropic, nematic, columnar, and crystalline phases appears, depending on the aspect 
ratio [l3|, [H| . Despite several previous studies (TH-Qij , the influence of electrostatic interactions on phase behaviour 
is still poorly understood and a matter of debate. Previous simulation studies of charged disks suffer from not taking 
into account the renormalised effective charge and the full orientational dependence of interactions. Furthermore, 
most of them are based on coarse-grained interaction site models [16-18], which is computationally costly. Here, such 
shortcomings are circumvented: our approach is based on a well controlled effective electrostatic potential between 
two disks, that depends on their relative distance and orientation. We shall neglect van der Waals interactions, size 
and charge polydispersity to focus on repulsive electrostatic interactions that we deem most relevant. We address two 
important questions: i) how do electrostatic interactions influence the isotropic-nematic (I/N) transition? ii) can the 
experimentally-observed arrested states occur within our model and what are their natures? 

We account for electrostatic interactions between charged disks, in an electrolyte medium with inverse Debye length 
ft, at the two-body level using Poisson-Boltzmann formalism [2(J, [21(. The far-field potential for two charged disks 
with diameter a, centres at a distance r apart, and with normals making angles 6i with the centre-to-centre line, 
reads: 
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FIG. 1. Characteristic anisotropy of electrostatic potential of charged disks in an electrolyte, a) Anisotropy 
function /(ft<7, 0) versus the angle between the disk normal u and the line of observation, for several values of ionic strength kg 
b) Two interacting parallel charged disks in stacked and coplanar configurations with equal centre to centre distance correspond 
to minimum and maximum repulsion, respectively. 



where Z e ff is the effective charge; it depends on the platelet bare charge and the ionic strength through k. For highly 
charged platelets as considered here, Z e ff saturates to a value Z^Xb/g = 0.5/^7 + 1.12, where Xb is the Bjerrum 

length (7.1 A in water at room temperature). The anisotropy function /(kg, 6) encodes the orientational dependence 
of the potential and can be expressed in terms of spheroidal wave functions (2l|. f(n<j,6) increases with angle 0, 
reaching a maximum in the disc plane (6 = it/ 2) (Fig. QJi). It furthermore is enhanced by the ionic strength: a 
smaller Debye length acts as a probe which reveals asphericity. This results in a strong asymmetry between coplanar 
and stacked configurations (Fig. [TJd), a notable difference with standard Heisenberg-like models. A competition 
between the excluded volume and electrostatic effects ensues, which is a generic feature for charged oblate spheroids. 
This interplay has important consequences on the phase behaviour and allows for the formation of new structures. 

We have investigated the structure and dynamics of hard disks interacting by the anisotropic Yukawa potential of 
Eq. [1] by Monte-Carlo simulations [22|. The degree of orientational order was characterised by the nematic order 
parameter S = ^(cos^)) = (!(3cos 2 ^ — 1)) where ijj is the angle of a platelet normal with the director n and 
the brackets mean averaging over all particles. Hence, — 1/2 < 5 < 1 with S = for the disordered isotropic phase. 
When < S < 1, the platelet normals point on average along the director n while in the more unconventional 
S < regime, the platelet normals are on average perpendicular to n (" antinematic" order). Upon varying density 
and ionic strength, we found eight distinct structures (Fig. [2j). Four of them have a vanishing S and are formed 
at low and moderate densities. Three of these structures, isotropic fluid, plastic crystal (structure with positional 
order but no orientational order), and "random stacks", are orientationally disordered while the fourth one is a novel 
structure that consists of intertwined nematic and anti- nematic layers (Fig. [2]). We coin it "nematic-antinematic 
intergrowth texture". The orientationally ordered phases with S > 0.4 which appear at higher densities include two 
liquid-crystalline phases, i.e. nematic phase of platelet stacks (sometimes called "columnar nematic") and hexagonal 
columnar phase, together with two crystalline phases with bcc-like and hexagonal symmetry. 

Before turning to a comparison with experiments, we present the phase diagram. It is essential to emphasise that 
our equilibration procedure for p* = pa 3 > 3 requires a careful and gradual increase of disk charge from zero to the 
desired value, to ensure attaining configurations close to equilibrium. More rapid quenches generally put the system 
out of equilibrium, and for high enough densities, we therefore expect discrepancies with the experimentally observed 
states, that have presumably not reached equilibrium. We now highlight important aspects of the phase diagram 
(Fig. [3]) in connection with the features of our pair potential. A significant hallmark is the non-monotonic behaviour 
of orientational disorder-order transition with kg which results from the opposing effects of the decreasing potential 
range and the increasing amplitude of the anisotropy function, upon increasing ionic strength; Z^f also increases 
with ionic strength but we verified that the observed trend is unchanged provided Z e ffXs/G > 2.3. 

In the limit of low ionic strengths (kg = 1), where the potential is rather isotropic and long-ranged, bcc-like 
structures are formed, as observed with charged spheres at low kg [23]. At low densities, disks are orientationally 
disordered (plastic crystal) and are reminiscent of Wigner crystals observed for low volume fractions of charged spheres 
[24j . Increasing the ionic strength, at kg = 2 where the amplitude of the anisotropy function is less than 1.2, the 
plastic crystal is replaced by an isotropic fluid at low and moderate densities. The crystal disappearance for a slightly 
larger kg and shorter range of potential highlights the geometrical frustrating effect of anisotropy. At kg = 2, in the 
high density regime, the long-range positional order is still preserved as crystals with hexagonal symmetry appear. 
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FIG. 2. Summary of found structures for charged disks. We recognise eight distinct configurations, four of which 
have a vanishing orientational order parameter S. These include isotropic fluid, plastic crystal, "random stacks" and a novel 
"intergrowth texture" that is composed of disks arranged in alternating nematic (red colour) and antinematic (blue colour) 
layers. The orient at ionally ordered structures with S > 0.4 encompass nematic and columnar hexagonal liquid-crystals in 
addition to two crystalline phases with bcc and hexagonal symmetry. Note that the disks are not displayed with their real size, 
but with 20 to 40 % their diameter. 



Interestingly, for large enough ionic strengths (4 < kg < 6) and at moderate densities, the new intergrowth texture 
appears and, for kg > 8, randomly oriented stacks of disks are observed. The intergrowth texture consists of sets of 
aligned disks (0.75 < S < 0.93) interspersed with layers exhibiting anti-nematic order (—0.45 < S < —0.3). Both 
types of layers share the same director. Here, the particles organise with relative T-shape configurations, intermediate 
between coplanar and stacked. At still higher densities, further increase of ionic strength (kg > 4), leads to weakening 
of the positional order. Hexagonal crystals are replaced by hexagonal columnar liquid-crystals with no positional 
correlations along the hexagonal axis. 

To better characterise the orient ationally disordered structures, we examined their translational and rotational 
dynamics by dynamic Monte-Carlo method [25j. We computed self-intermediate scattering functions F s (q,t) = 
(exp[iq - (f(t) — f(0))]) and time orientational correlation functions (P2(u(t) ■ u(0))). In the isotropic fluid, both 
translational and rotational correlation functions decay on a rather fast time scale (Fig. 2]), whereas, in the random 
stacks phase and intergrowth texture, the orientational time correlation functions show very little decay within the 
longest time-scale probed, in contrast to their translational counterparts. In random stacks, increasing p* from 2 to 
3 has a spectacular effect on orientational relaxation (Fig. [4}d), which becomes completely hindered; the dynamics of 
this structure that lacks long-range orientational order, is reminiscent of that in orientational glasses [26|. This state 
results from the strongly anisotropic and short-ranged character of the potential for sufficiently large kg, where the 
locally favoured stacking leads to frustration and hinders the onset of orientational ordering. On the other hand, the 
plastic crystals exhibit slow dynamics for translational coordinates and fast decay for orientational ones. 

We now compare our simulations against the most extensively studied charged colloidal platelets: Gibbsite, Bei- 
dellite, and Laponite systems. Gibbsite and Beidellite show an I/N transition for a wide range of ionic strengths, 
and platelet stacking is often observed in the nematic phase. Moreover, Gibbsite suspensions also display a columnar 
hexagonal liquid-crystalline phase. These features are well captured in our simulations. Fig. [5] shows the phase dia- 
grams of Gibbsite [4] and Beidellite [3], drawn versus reduced density p* and kg for quantitative comparison. Overall, 
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FIG. 3. Phase diagram of repulsive charged disks as a function of dimensionless density p* = pa and ionic strength 
kg. The vertical dashed line shows the density beyond which the nematic phase for hard disks (thus without electrostatic 
interactions) appears. The solid line separates the orientationally ordered (right hand side) and disordered phases (left hand 
side). The dashed line is a guide to the eyes to show, on the right hand side, the regions where the structures with slow 
orient at ional dynamics appear. 




FIG. 4. Dynamics in isotropic fluid and random stacks structure, a) Self-intermediate scattering function F s (q,t) 
computed at qa — 7 and b) orientational time correlation function (P2{u(t) • ^(0))) for an anisotropic fluid kg — 20, p* = 2 and 
random stacks kg — 20, p* = 3. The time-scales are expressed in terms of characteristic time of free translational tb = cr 2 / (6Dq) 
and rotational tq = 1/(2 Do) diffusions, respectively. 



the agreement between our simulations and the experimental phase diagrams for Gibbsite [4j and Beidellite [3] (shown 
in Fig. [5j) is quite remarkable. In both Gibbsite and Beidellite, for 10 < kg < 20, a nematic phase of stacked platelets 
is found, although at densities slightly lower than the simulated ones. This small discrepancy could be due to platelet 
polydispersity that usually widens the isotropic-nematic coexistence region. Furthermore, for Gibbsite particles at 
lower ionic strengths, a direct transition from isotropic to columnar hexagonal phases is observed, which is again 
similar to our results at kg = 8. Upon further increase of density, a nematic gel is experimentally observed whereas 
we obtained a columnar hexagonal phase. At such high densities, we already anticipated that the suspensions cannot 
easily reach equilibrium, since the simulations require a slow and careful procedure to relax. Furthermore, the nematic 
gel of Gibbsite is a jammed state of oriented stacks (columns) that evolves towards a columnar hexagonal structure 
with time p7j . 

In contrast to Gibbsite and Beidellite, the I/N transition in Laponite suspensions is hindered by the occurrence 
of various kinetically trapped states at very low volume fractions (<j) < 0.02) and formation of these states is shifted 
towards lower densities upon increasing ionic strength The range of densities and ionic strengths investigated in 
experiments corresponds to p* < 1 and 1 < kg < 20. At such low densities, we observe a solid-like state (plastic 
crystal) only at small kg = 1. This state could presumably only be glassy if polydispersity were considered (28) . 
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FIG. 5. Comparison to experiments. Experimental phase diagrams of a) Gibbsite [|| and b) Beidellite .3]. The phases 
corresponding to each symbol are indicated in the legend. Here, I refers to isotropic, N to Nematic stacks phase and NG to 
nematic Gel. 

Upon further increase of ionic strength, we observe an ergodic repulsive isotropic liquid in the low-density region of 
our phase diagram. Interestingly, at somewhat higher densities and for kg > 10, we find signatures of an orientational 
glass that is formed at lower densities upon increase of ionic strength, similar to the trend observed for Laponite. 
Our model system only based on repulsive interactions cannot fully explain the phase behaviour of Laponite for 
larger ionic strengths; these discrepancies point to the possible relevance of attractive interactions [29]. Nevertheless, 
our approach hints at the existence of a solid-like state at low ionic strengths, in agreement with the experimental 
observation of a Wigner glass [30| . 

In conclusion, numerical simulations of charged platelets, interacting with an effective potential obtained from 
Poisson-Boltzmann theory, show a rich phase diagram. We thereby capture the generic features of the experimental 
systems. The required procedure of slow charge increase, at moderate and high densities, echoes the difficulty for 
experimental systems to reach thermodynamic equilibrium. Better agreement with experiments could be achieved 
by considering polydispersity, platelet flexibility, and interactions specific to each system such as van der Waals 
attractions. The features brought to the fore here, though, can be considered as the non-specific effects pertaining to 
charged colloidal platelets. In particular, the robust form of our screened Coulomb potential leads to the formation of 
orientational glasses (random stacks and novel intergrowth texture), which calls for further experimental investigations 
in the isotropic region of the phase diagrams close to the I/N transition. 

I. METHODS 

We carried out canonical (NVT) Monte-Carlo simulations 22] on a system of N thin hard disks of diameter a in a cubic 
simulation box of side L with periodic boundary conditions using the anisotropic two-body potential of Eq. [TJ The simulations 
were performed at T = 300 K for a wide range of densities p* = Ng 3 /L 3 = 0.1 — 8 and ionic strengths kg — 1 — 20. The diameter 
g/Xb was taken to be 43.1 in the simulations. The number of particles varied between 500 < N < 10976, although most of the 
simulations were performed for N = 1024 or N = 1372. For ionic strengths kg < 2, we adopted an Ewald-like scheme to take 
into account the relatively long range of the potential [3~0 |. Due to the rather large values of the effective charge, we encountered 
some convergence difficulties for equilibration of large densities, p* > 3 as the obtained final configuration depended on the 
initial one. Therefore, we adopted the following procedure: we started from an initial equilibrium configuration of (uncharged) 
disks at the desired density; then increased the charge gradually to the final value Z^^kg) allowing the system to equilibrate 
at each step. The required procedure is akin to simulated annealing used in spin glasses. 
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